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Abstract. In this note, we discuss dilation-theoretic matrix parametrizations of 
contractions and positive matrices. These parametrizations are then applied to some 
problems in quantum information theory. First we establish some properties of posi- 
tive maps, or entanglement witnesses. Two further applications, concerning concrete 
dilations of completely positive maps, in particular quantum operations, are given. 



I. Introduction 

It is well-known that positive operator-matrices, or more generally positive kernels, 
can be parametrized by contractions 5J. In this paper, we show that analogous results 
can be obtained for matrix contractions. A common feature of these parametrizations 
is that, while the explicit expressions may seem intricate, due to their combinatorial 
nature they can be easily understood by inspecting the so-called lattice diagrams. 
These diagrams will be used repeatedly to illustrate accompanying results. 

The organization of this note is as follows. The structure of row and column con- 
tractions are discussed first. They already possess an elegant combinatorial structure 
and play a central role in our parametrizations. Next we consider matrix contractions. 
The 2x2 matrix contractions were already characterized in |2|. Here we extend the 
description to matrices of arbitrary size and point out the combinatorial aspect of this 
parametrization. Then the special case of unitary matrices is examined. We also review 
the parametrization of positive matrices. While the definitive treatment of the positive 
case is [2], our discussion differs slightly in some minor technical details. Turning to 
applications, we first obtain some properties of positive maps. The structure of positive 
maps and contractions is applied to show that general positive maps are more than 
merely positive when restricted to certain subsets of positive matrices. Results of this 
type were obtained in 0] and the parametrization of positive matrices allows one to 
explore their extensions in a non-ad hoc way. By the correspondence between positive 
maps and entanglement witnesses [7j, we thus show that certain families of bipartite 
mixed states are separable. The last two applications concerns the unitary dilation of 
completely positive maps on matrix algebras. While a celebrated result by Stinespring 
[TT] shows that such dilations always exist, the paramatrization of contractions allows 
one to give a concrete constructive procedure for such dilations. 



The results on positive maps are first stated in operator-theoretic terms before being 
placed in physical context. The last two applications are phrased more directly in 
the language of quantum information theory. For general background in quantum 
information, we refer the reader to and 

Although we only consider matrices of finite size, all parametrizations described in 
this paper can be extended to (semi-)infinite matrices, where convergence is given by 
an appropriate operator topology. 

2. Row and Column Contractions 

The facts outlined below can be found in [Sj, in which the structure of such con- 
tractions plays a crucial role in obtaining a parametrization for positive kernels, the 
Schur-Constantinescu parameters. Our presentation here is different in that the notion 
of defect spaces is dispensed with. In [5 , the defect spaces are used in forcing the 
uniqueness of certain operator angles. Rather, we identify the unique positive square 
roots explicitly via partial isometries. 

In the following, H and Hi denote Hilbert spaces, and £(7ii,7i 2 ) the bounded op- 
erators from Hi to H2- 

Lemma 1. Let X G C(Hi,H2) and Y G £(Hi,H2) be bounded operators between 
Hilbert spaces. X*X < Y*Y iff there exist a contraction T : H2 — > H2 such that 
TX = Y. 

Proof: Suppose X*X < Y*Y. Define T' : by T'Xh = Yh on RanX, the range of X. 
Extending by continuity to the closure of RanX and then by to the orthogonal com- 
pliment gives a contraction T. The expression X*X < Y*Y implies KerY C KerX. 
So T is a contraction satisfying FX = Y. The converse is also straightforward. □ 

This leads to the following well-known fact regarding the freedom in the square-roots 
of a bounded operator. 

Lemma 2. If X*X = Y*Y , then there exist a partial isometry V such that VX = Y. 
Equivalently, X*V* = Y* . 

Proof: It is clear that we can take V to be the T defined above. □ 

The partial isometry V is unique if the condition KerV C (RanL) 1 - is imposed. 
This uniqueness condition will be assumed throughout. The operator V is unitary if 
both RanX and RanY are dense. When Hi = H2 = H3 = H is finite dimensional, 
V can aslo be assumed unitary. Given a positve operator A G £(H) with unique 
positive square root A2, every L such that L*L = A is related to A 2 by A2 = VL, or 
A\ = L*V*. 

These two lemmas will be used to exhibit the structure of row contractions. Before 

doing so, we first introduce a bit of terminology. For a contraction T G C(Hi,H2), the 
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positive operator (I — T*T)i G C{7i\) is called the defect operator of T and denoted 
by D T . Similarly, D T * = (I - TT*)^ G C(H 2 ) is the defect operator of T*. When 
Hi = H 2 , by the continuous functional calculus, 



implies TD T 



T{I - T*T) = (I - TT*)T 
Dt*T. This relation shows that the operator 



J(T) 



T Dt* 
D T -T* 



, called the Julia operator of the contraction T, is unitary. The Julia operator can be 
viewed as the building blocks of our parametrizations. It is represented by figure 1. 




Figure 1 . The Julia operator J(T) 

Next we describe row contractions of length two. 

Proposition 1. Let T = [7\ T 2 ] G C(Hi ®H 2 ,H), then ||T|| < 1 iff there exist 
contractions Y\ G £(7ii,7i) and T 2 G C(H 2 ,H) such that T = [Ti ^^2]- 

Proof: (=>•) We can take r : to be T x . \\T\\ < 1 implies 

j _ TT * _ j _ rir * _ t 2 t* > 

i.e. Dp* > T 2 T 2 *. By lemma, AD r * = T 2 *. choosing T 2 = A* finishes the argument. 
(<=) Direct computation. □ 

The defect operators for the contraction T = [Ti T 2 ] can be directly calculated: 



D 2 T = 

Invoking lemma, we have 



D 7 



D Tl 

-r*r x d Ti 

for some unique partial isometry V . Similarly, 



V 



D 2 T , = (D n D n D n D rt ). 
The general structure of row contractions is described by: 

Theorem 1. Let T =[T 1 T 2 • • •] G C{®U h H), then \\T\\ < 1 iff there exist contrac- 
tions Ti G C(Ti,i,H) such that 

T=[T 1 ,D ri T 2 ,D r *D r *T 3r -- , D r * ■ ■ ■ D r * n _r n \. 
Furthermore, the defect operators are of the form D\ = 



-r*r x 



o 

D T2 



T* -Dr* ' ■ ■ Dr*Fi — r*-Dr* ' ' ' Dr*^2 








D T2 



-r*D r * 

J- 1 2 

— r^-Dr* 



and 



-D^» — Dr* ■ ■ ■ -Dr* -Dr* ■ ■ ■ -Dr* . 

1 1 1 1 n 1 n 1 1 



Proof: The argument is by induction. The length n = 2 case was shown above. Now 
suppose the claim holds for length n — 1. For a row contraction T of length n, put 



T=[S,T n ] 



where by inductive hypothesis, 



s = [r\, D r *r 2 , Dr-D^rs, • • • , D rr • • r n _!] 

According to lemma, there exists a contraction A such that 



T n = D S *A. 



But 



D 



D q* — Dr* • • • Dr* V 

for some partial isometry V. Choosing T n = VA shows T is of the desired form. 
Applying the defect operator identity proves the remaining proves the remaining claims. 
□ 

The combinatorial content of the theorem can be depicted pictorially. Figure 2 be- 
low shows the parametrization of length 3 row contractions. The downward arrows 
indicate input ports and the upward arrows output ports. For example, for a matrix 



(Tij), each path from input 3 to output 1 gives rise to a term in the expression of the 
T 13 entry (in this particular case, T 13 is the third entry of a row contraction T). 




Figure 2. Structure of row contractions of length 3 

Figures 3 and 4 describe the parametrizations of the defect operators (more pre- 
cisely that of the natural square roots, or Cholesky factors). In the special case that 
TT* = I n , i.e. T is a surjective partial isometry, the contraction T n is in fact a partial 
isometry. 






FIGURE 3. Structure of D T where T is a row contraction of length 3 

The structure of column contractions can be exhibited in similar fashion. For com- 
pleteness, we state the corresponding result for column contractions. 

Theorem 2. An operator 



T 

L ± n 



is a contraction if and only if 





^1 




Figure 4. Structure of where T is a row contractions of length 3 



Ti = Ti and T k = TkDr^ ■ ■ ■ D Tl , 
where IVs are contractions. Moreoever, the defect operators ofT and T* take the form 



and 



D 2 T = D Tl ...D Tn D Tn ...D Tl , 



Drp* 



-r 2 r? 



o 

1 o 








D r2 Tt -T n D T „ 



D 



p* 








-r\r* 





-ri^r, 

-r 2J Dr 3 



One can be easily convinced that there are lattice diagrams corresponding to the 
above description. If a column contraction T is such that T*T = I n (i.e. T is an 
isometry from Ti to T n is a partial isometry. 



Let T : Tii ® Ti 2 —> KLi 



3. Matrix contractions 

^2, 



T 



A B 
C D 



be a contraction. Then L4 B\ = P^T, where P^ denotes projection onto /Ci, is nec- 
essarily a row contraction, therefore of the form [ri -D r *T 2 ] where Tj are contractions. 
Similarly, we have 



" A ' 






C 




_ r 3J D ri _ 



Dr 



So we have 



T 



T 3 D ri D 



We will show that the entry D, unspecified so far, can also be parametrized by con- 
tractions. To this end, view T as a row contraction T = [Si S2] with 



Si 



T 3 D ri 



and So 



D r *r 2 

D 



Therefore S2 = D$*A for some (column) contraction A. S\ is a column contraction 
and direct calculation gives 



Thus 



So 



Dri 



D vl 



-v x v% 



D 



r 3 



Ai 
A 2 LV 



Comparing entries and invoking the uniqueness condition shows that Ai = T 2 . Rename 
A2 as T4 and we have 



D = -r 3 rjr 2 + D r ,r 4 D r3 . 

The above can be summarized by j2]: 
Theorem 3. Every 2x2 contraction T : 7ii © 7i 2 —> K,\ © /C 2 of the form 



T 



where T,- are contractions. 



Ti L> r *r 2 

t 3 d Fi -r 3 rir 2 + D T3 r A D T2 



The general structure of a n x m matrix contraction can be obtained in a similar 
way. 

Theorem 4. Lei 



11 



L T nl 



T 



\ Hi > ©1 K-i 



be a contraction. Then the column contraction 

7 



L ^nl J 



is of the form 



r 2 £> ri 



iV*^ • • • d Fi 



where Tj are contractions. 



For 1 < k < m, the column 



ik 



T 

. ± nk . 



is obtained inductively by 



ik 



T 

. ± nk . 



Dc^Dc* ■ ■ ■ Dci__C, 



where Ci is the column contraction parametrized by r n (j_i) + i...r ni with C\ being as 
specified above. 

The proof is immediate and omitted. The combinatorial structure of matrix contrac- 
tions can also be nicely described via lattice paths. For example, the lattice diagram 
for the 2x2 case is figure 5. 




Figure 5. A 2 x 2 matrix contraction 

The defect operators of matrix contractions can also be calculated. Due to the "two- 
layer" nature of its parametrization, the explicit formulae may seem complicated. It 
is helpful to first look at the lattice diagrams. From inspecting the above figure, we 
anticipate that the Cholesky factor of D T and D T * to have the corresponding pictures 
given by figures 6 and 7 respectively. 

In other words, one should have 



Dn 



D F3 D ri 

-r^Drs - £r 2 r:r 3 d T2 d Ta 



D Tz D Tl -D T . A Y\Y 2 - YIY A D T2 
D r4 D r . 2 



and 
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Figure 6. The natural Cholesky factor of Dt where T is a 2 x 2 matrix contraction 
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1 
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Figure 7. The Cholesky factor for D T * 



D n D T% 

This can be confirmed by a straightforward but perhaps tedious calculation, which we 
shall not bore the reader with. The defect operators for matrix contraction of any finite 
size can obtained in similar fashion. 



1 1 1 2 







-r 2 r* 4 D r * - D r *r*r* D r *D n 



Unitary Matrices The unitary operators are the extreme points of contractions, thus a 
special case. If 



: H x © H 2 -> /Ci © K, 2 



T= ' Ti D r .r 2 

[ t 3 d Ti -r 3 rjr 2 + d T3 v a d V2 

is unitary, then [I^ D r *T 2 ] is a partial isometry therefore so is T 2 ; same goes for T 3 . 
When all spaces are finite dimensional and, in the expression 



T = 



A B 
C D 
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we have B and C being square matrices, T 2 and 1^ are unitary and the description of 
T becomes very simple: 



T 



Ti D T *r 2 

V 3 D ri -T 3 TIT 2 

In other words, all unitary matrice are related to the Julia operator via 



T 



' 1 










r 3 . 





ri 

An 



1 1 





" 1 










r 2 . 



This will be applied in the sequel in calculating dilations of completely positive maps/ quantum 
operations. 

4. Positive Matrices 

Similar dilation-theoretic techniques can be applied to positive matrices. As stated in 
the introduction, we outline basic results for completeness. See jS] for a comprehensive 
discussion. As for contractions, one can start by examining the 2x2 case then apply 
induction. Let 



.4 



G £{Hi@H 2 ) 



L\L\ A\2 
-^12 L 2 L 2 

be a positive semidefmite operator matrix whose entries are bounded operators, that 
is 



h 2 



,A 



h 2 



for all hi e TCi and h 2 ETC 2 . 



Theorem 5. There exists an unique contraction F G C(Ti. 2 ,Ti.i) such that A i2 = 
L\YL 2 . 

Proof: : Assume for the moment that both An and A 22 are invertible. Then a 
Frobenius- Schur identity holds: 



A 



L\Li A12 
Air. LXLo 



I 

A\ 2 {L\Li)- 1 I 



L\Li 
L\L 2 -A\ 2 {L\LiY l Ai 2 



l 12 -^2-^2 

It follows that A is positive if and only if its Schur complement 



I {L\Li)~ 
I 



L^-Al^LlLi^Au > 



i.e. 



L* 2 L 2 > A\ 2 {Li)-\L*i)- l Ai 2 . 
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By lemma 1, there exist a contraction T such that TL 2 = (L 1 l )*Ai 2) where we impose 
the condition required for uniqueness. Thus Ai 2 = L\TL 2 . 

For the general case where An = L\L\ and A 22 = L 2 L 2 need not be invertible, 
consider the sequences {a*a n = A u + ^} and {/?*/3 n = A 22 + -}■ By the spectral 
mapping theorem for self adjoint operators, a* n a n and /3*/3 n are invertible for all n. 
Therefore there exist contractions {r„} with A\ 2 = a^T n /3 n . Since the unit ball in 
£(H 2 , Hi) is compact in the weak operator topology, T n converges to some contraction 
T weakly We can compute directly, for all hi 6 Hi and h 2 e H 2 , 

(hi,LlTL 2 h 2 ) Hl = lim(hu L* 1 T n L 2 h 2 )Hi = ^{h,a^T n p n h 2 ) Hl = (h, A 12 h 2 ) Hl . 

n n 

This proves the claim. □ 



This can be generalized to positive operator matrices of arbitrary size in the obvious 
way. We present the finite case as an algorithm 

Algorithm 1. [Sj Let A = [A^]^ E £(©" = i7"^) be positive. The Schur-Constantinescu, 

or SC, parametrization of A can be calculated recursively as follows: 



An—l^n—l A n —i^ n 

A , A 

■ rl n,n—l ^Mi,n 

the 2x2 case. 



i) 



is positive and can be parametrized according to 



ii) For 1 < k < n — 2, the SC parametrization of 



is calculated by first considering 



A 



fc+i,fc+i 



A 



n,k+l 



Ak,k 


Ak,k+l 


Ak,n 


Ak+i,k 


Ak+i,k+i 


Ak+l,n 


An.k 


A n ,k+i 


A 


Ak+l,n 








— T* 


-l-^fc+l,fc+l; 


A 







where L k+ i is the Cholesky factor calculated in the previous step then put 



[Ak,k+i ■ ■ ■ A k:n } — Al k R k Lk+i 
with R k being the corresponding row contraction. 

The lattice diagram for a 4 x 4 positive mattrix is given below. 



Tensor product of positive matrices If M = (m^) 6 Q nxn j s a positive matrix with 
complex entries and A = B*B 6 £(H) where H is a Hilbert space. Suppose M is SC- 
parametrized by {7^} C C, i — 1, - • • , |m(m — 1) with | -y^ | < 1. The parametrization 
of the positive matrix 



M ® A = (m i3 L*L) e C nxn g> C(H) 
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Figure 8. Lattice structure for 4x4 positive matrices 



can be described easily jHj- Namely, let Tj = jJ-h and Lj = y/muB\ it is clear that 
they parametrize M (g) A in the sense of Schur-Constantinescu. 



Matrices given by a strict inequality. The natural square roots given by the SC- 
parametrization are upper(or lower)- triangular, i.e. they are Cholesky factors. Owing 
to this fact, if A* A > B*B, and the SC parameters of A* A are known, B is readily 
described. Take for instance the 2x2 case. B = FA, where T is a contraction. Let 

ri D r *r 2 
t 3 d Fi -r 3 r*r 2 + d T3 t 4 d T2 J ' 

then TA corresponds to the following figure: 



A 







AL 22 



and T 




5. Applications 



Due the the ubiquity of positive matrices, the Schur-Constantinescu parametrization 
of positive matrices has numerous applications [5]. More recently, it has been applied 
in the context of quantum information theory. For example, it was used to parametrize 
completely positive maps (quantum channels) in A cylinder-like condition, called 
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the Bloch cylinder, was obtained for positive matrices of trace 1 (quantum states) of 
any finite dimension. This provides an alternative to the well-known Bloch sphere. 
In |T2] it was applied to show that every positive map is completely positive to a 
certain extent, thus establishing the separability of certain families of quantum states 
in arbitrary finite dimensions. In a similar vein, in this section we obtain more results 
in this direction, in a sense extending what was found in Also, we consider two 
further applications that are matrix completion problems in disguise and can be solved 
by utilizing parametrization of matrix contractions. 

5.1. Positive Maps. In this section, the structure of contractions and positive matri- 
ces are applied to extend properties of positive maps. 

Definition 1. Let 7i and JC be Hilbert spaces. A linear map $ : C{7i) — > C(IC) is 
said to be positive if it preserves the cone of positive elements, i.e. A > implies 
> 0. Let C nxn denote the n x n matrices of complex numbers and I n the identity 
map on C nxn , then a map $ is said to be n-positive if the induced map 

I d ®$: C nxn ® C(H) -> C nxn ® C{K) 
is positive, and $ is completely positive, or CP, if it is n-positive for all n. 

We state the following result without proof [TU] . 
Theorem 6. (Russo-Dye) Let <3> be a positive map between unital C*-algebras, then 

\\n < \w)\i 

In particular, if $ is unital and V a contraction, then 

$(r*)$(r) = ||$(r)|| 2 < ft < i. 

Similarly, <3>(T)$(r*) =< /. Making use of this, one has 0: 

Theorem 7. Let $ : C(7i) — > C{K) be a postive map, then for all positive 



T S 
S* T 



eC 2x2 ®£(H), we have (I 2 ®$)(A) 



$(5*) $(T) 



> 0. 



We recast the proof so that the role played by contractions is made more apparent. 
Proof: Assume for the moment that T _1 exists and is invertible, therefore so is 
$(T). According to theorem 5, S — TaTTa for some contraction T. It is equivalent to 
show that the Schur complement 



$(T) - $(5*)$(T)- 1 $(5) > 

i.e. 



I > <$>{T)-^§(T^T*T^)§(Ty^ ■ $(T)"^$(T^rT^)<l>(T)"i 
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This suggests that we define \I/ : C{7i) — > £(/C) by 



V(A) = $(T)-5$(T5AT5)$(T)-i 
\1/ is an unital positive map. By Russo-Dye, 

/ < # (r)tf(r*) = $(t)-^(t^t*t^)$(t)-^ ■ $(r) - 5$(T5rT5)$(r)~5, 

which is what we want. 

If T is not invertible, consider the sequence {T n = T + -}. T n tends to T uniformly 
and positive maps are bounded. So 

$(5)=lim^(T n )3-A n .li m $(T n )3 

n n 

for contractions A n . Let A G £(/C) be a weak operatorial limit of A„, then 

= $(limT n )5 . A . $(limT n )5 = $(T)U$(T)i 

n n 

So the claim holds. 

If, in addition, $(/) is not invertible, take a positive functional (p with 0(7) = 1. 
Define $ n (A) = $(A) + -0(A). We have $ n — > $ in the operator norm of linear maps, 
and 



$(S) = lim$ n (S) = lim$„(T)5 . A^ ■ lim$ n (T)3 



The same weak limit argument shows that $(<5) = $(T) 2 A'<I>(T) 2 for some contraction 
A'. This proves the theorem. □ 

In other words, any positive map is 2-positive on the 2x2 Toeplitz matrices. Now we 
extend this to a sub-family of 3x 3 positive matrices. Recall that an operator A G C(TC) 
is said to be subnormal if it is the compression of a 2 x 2 normal upper-triangular N, 
i.e. if there exist some Hilbert space K, and a normal N G C(K) such that N is of the 
form 

" A B 
C 



N 



The following fact, which we state without proof, will be used |3] : 
Lemma 3. For any unital positive map $ : C(7i) — > C{K) and any normal A G C(7i), 



<5>(A*A) > $(A*)$(A) and <$>(A*A) > $(A)$(A*). 

What is known as Kadison's inequality will also be needed: for every unital positive 
map $ and every self-adjoint S, $(5' 2 ) > "^(S*) 2 What we will show that is essen- 
tially every positive map is 3— positive in a certain limited sense. We first notice that 
subnormal contractions enjoy a property stronger than that prescribed by Russo-Dye. 
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Lemma 4. If $ is a unital positive map and T a subnormal contraction, then 



I- $(rT) - $(D r «) 2 > and I — $(TT) — $(L> r ) 2 > 0. 
Proof: We directly compute 

/ - $(rT) - $(D r ,) 2 > / - $(r*r) - $(£>jl) 
= / - $(r*r) - $(/) + $(r*r) = $(rr) - $(r*r), 

which is positive, by the preceding lemma. The second inequality is similar. □ 

Theorem 8. i) Consider positive matrices in {A^) G C 3x3 <E> that are SC- 

parametrized in the following way: for i = 1,2,3, choose An = T > 0. Of the three 
contractions, choose I\ 2 to be subnormal, r 23 = 0, and I\ 3 = /. In other words, we 
consider 3x3 positive matrices of the form 



(A 3 ) = 



T T2TT2 TzD r *T2 
T 
T5D r *T5 T 



Then 



(j 3 <g>$)(4i)>o, 

for any positive map $ acting on £(H). 



(4,0 = 



We do not completely recover Choi's result by considering the 2x2 leading minor 
in part i). The requirement that T be subnormal is particular to the 3x3 case, due to 
the presence of D^* in the parametrization. 

Proof: i) Assume for the moment that 3>(I) and T are invertible. Invoking again 
theorem 5 on the 2x2 case, it is equivalent to show that 



is subnormal, T 


12 — 


0, r 13 = /, 


T 










T 








\ T 



$(T) > [ $(T*TT*) $(T^D r .T2) ] 
The right hand side is 



^(T)- 1 
$(T) _1 



$(T2r*T2) 
$(T5£) r *T5) 



$(T^rT5)$(T) _1 $(T^r*T^) + $(T*D r *T*)$(T)- 1 $(T*D r .T*). 
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Again we define an unital positive map \l/ by 

#(£) = $(T) _ 3$(T^BT^)$(T)-^. 
By the subnormality of T and lemma 4, 

/ > *(r)*(r*) + ^(D r *) 2 

whicn proves the claim. The general case can be shown using argument similar in that 
of theorem 7. 

ii) The argument is analogous to i) and omitted. □ 

The following result is of similar nature. The special 2x2 case, proven in says 
that any positive map is positive on 2 x 2 Hankel matrices. 

Theorem 9. If A is a positive matrix of the form 





' T 


• 


• 


Si ' 




' T S 1 ■■■ 




S m -i 









T • 


• 


s 2 




Si R 









A = 










or 


: "-. 






e C mxm ®£{H) 







• 


• T 








R 









Si 


S 2 • 


■ 


R 




S m -i o 





R 





where m is arbitrary, then 



for any positive map $. 



(I m ®$)(A) >0 



Proof: Suppose A is of the first form in the claim. The argument uses only the structure 
of 2 x 2 matrices and thus we consider first the case m — 2. By virtue of previous 
arguments, it can be assumed without loss of generality that T and <&(/) are both 
invertible. To show 



$(T) $(S) 



> 



, it suffices to obtain $(ST 1 S') > $(S')$(T) because the Schur complement 

R — ST^S is positive. This is equivalent to 

$(T)-5$( ( 5T- 1 5)$(T)-5 > $(T)-s$(S)Q(T)-* ■ $(T)-?Q(S)$(T)-s, 
which suggests we define a unital positive map by 



V(B) = $(T)~2$(5)$(T)" 
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Invoking Kadison's inequality then proves the 2x2 case. For m > 2, it is enough to 
show that 



H[ Si 



T -l 







T -l 



S m -1 



i.e. 



m—1 



m—1 
i=l 



m—1 



i=i 



i=l 



But A is positive only if the principal 2x2 minors 

T Si 
Si R 

are positive. Thus the desired inequality holds by linearity of $ and the 2x2 argument. 
□ 

Remark What we have show above is that a positive map is 3-positive and CP (in 
the case of theorems 8 and 9 respectively) to a certain extent. Results in the similar 
vein were obtained in ^2] that are also applications of SC parameters. Namely positive 
maps were shown to be CP if restricted to certain families, of arbitrary finite size, which 
can be SC-parametrized by two real parameters. In that approach , Choi's result on 
2x2 Hankel matrices also become a special case. For comparison, a result from 
for matrices with entries in C 3x3 is stated below. 

Theorem 10. Let S be the linear span of 

a a b 

{A 



a,b,be C} c C 3x3 
and where m G N be arbitary. Then for all positive p G C mxm ® S, 



a a b 
b b c 



(/ w ®$)0o) > 

where $ is any positive map acting on C 3x3 . Furthermore, the claim holds if A is 
replaced by 



a 


c 


a 




a 


c 


c 


c 


b 


c 


or 


c 


b 


b 


a 


c 


a 




c 


b 


b 
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Separable Quantum States In physical language, trace-class positive matrices are un- 
normalized mixed states pQ. 

Definition 2. ^3] Let the state space of a bipartite quantum system be the tensor 
product H — Hi ®H2, where Hi are Hilbert spaces. A state a 6 £(H) is separable if 

it lies in the closure, in the trace norm, of states of the form 

k 

where p\ are states in Hj . 

The membership problem for separable states is sometimes called the separability 
problem. The following theorem establishes the correspondence between the classifica- 
tion of positive maps and the membership problem for separable states[7j: 

Theorem 11. If a mixed state a G C(Ha) ®C{Hb) is such that for every positive map 
$ from £(Hb) to C(TCa), the operator {I a ® < ^ ) )(o") is positive, then a is separable. 

The above result is of geometric nature and a consquence of the hyperplane-separation 
variant of Hahn-Banach. Thus in quantum information theory, positive but not CP 
maps are called entanglement witnesses, for they detect the entanglement of some state. 
If a family S of positive matrices is such that any positive map behaves as a CP map 
when restricted to S, then S must consist of separable states, due to lack to entangle- 
ment witnesses. Thus what was shown in the previous discussion translate to that all 
positive matrices of the forms specified in theorems 8 and 9 are separable states. In 
particular, 2 x m block- Toeplitz and block-Hankel states are separable, where m need 
not be finite. 

5.2. POVM's. We first give a few relevant definitions and basic results; the reader is 
referred to PQ for more background information. In the von Neumann measurement 
scheme, the effects of a quantum measurement are assumed to satisfiy the projective 
hypothesis, i.e. they are self-adjoint projections and form the resolution of identity 
of a self adjoint operator. A resolution of the identity is sometimes called projection- 
valued measure, or PVM. A more general formulation of measurement replaces these 
projections by positive operators: 

Definition 3. ^Uj Let X be a compact Hausdorff space, B be the Borel a— algebra on 
X , and H a Hilbert space. A positive operator-valued measure, or POVM is a 
map E : B — > C(7i) such that E(B) > for all B, and E is countably additive in the 
weak topology on Ti, i.e. for any pairwise disjoint collection {Bi}i>i C B, 

(E(U i B l )x,y) = J2(E(B l )x,y) 

% 

for all x,y 
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If E(B) is self adjoint for all B and E{B X n B 2 ) = E(B 1 )E(B 2 ), then each E(B) is 
a self adjoint projection. Thus E is a PVM and we recover the von Neumann Scheme. 
A natural question is whether a POVM can be "lifted" to a larger space where it is a 
PVM. The general answer to this dilation-theoretic question is given by: 

Theorem 12. (Naimark) Let X be a compact H aus dor ff space. Suppose E is a POVM 
on the o-algebra generated by the Borel sets of X and E takes values in C(7i). There 
exist a Hilbert space K, containing 71 as a subspace and a PVM F on X with values in 
£(/C) such that 



E(B) = P n F(B)P H 

Jor all Borel B C X. 

The proof follows from the fact that C (X) is a commutative C*-algebra and therefore 
the induced map 

/ G C{X) / f(x)dE(x) 
Jx 

is completely positive, rather than merely positive. Stinespring's theorem on CP maps 
jTTj then says $ can be dilated to a homomorphism The PVM corresponding to $' 
is the desired F. For a complete proof, see JU]- Stinespring's theorem is a generaliza- 
tion of the Gelfand-Naimark representation thereom of postive functionals. 



In quantum information theory, Of particular interest is the case when X is finite, 
with the discrete topology. In that case, one would like and can indeed find solutions 
of more concrete nature. Let X = {1,2,- ■• , n}. Without losing generality, we can 
consider only POVM's whose elements are rank-1 projections that may not be mutually 
orthogonal. Suppose a POVM on X is given by E(i) = ViV*, % = l,---,n with 
Yli v i v i = Im where m < n and I m is the identity in C mxm . In other words, 



M — [vi ■■■ v n ] G C mxn 

is an isometry, i.e. M*M = J. We want to specify a PVM F taking value in some 
C(K) whose restriction to C m is E. This is a trivial completion problem: given a 

(rectangular) isometry M, find a suitable rectangular N s.t. 1 



N 



is unitary. It is 



an elementary fact from linear algebra that such an N can always be found. 



However, in physical considerations, a suitable POVM is often obtained by coupling 
to the original system an ancilla. This amounts to finding appropriate operators A, B, 
and C such that 



U 



M A 
B C 
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G C 



kxk 



is unitary. The columns of U, {u±, ■ ■ ■ ,Uk}, then gives a PVM on C k = C m © C h 
with the desired properties. 



The completion of 2 x 2 unitary matrices offers an immediate solution to this problem. 
Since M is an isometry, it is trivially a contraction. For example, consider the special 
case of Julia operator. So the PVM are the projections onto the column vectors of 



J(M) 



M D M * 
D m -M* 



M*M = I means D M = and similarly D M * = (I - MM*] 
I — MM* is a projection. Thus 



[I - MM*) because 



J{M) 



M I 




MM* 
-M* 



and a suitable PVM F is obtained without doing any calculation whatsoever. If M is 
mx n with m < n, the Hilbert space K, is of dimension m + n. From the discussion on 
unitary operators, we see that the freedom in obtained a suitable PVM is described by 



J'(M) = 

where Ui are unitary matrices of suitable size. An effect of the POVM can be obtained 
by first performing the corresponding yes-no measurement from the dilated PVM then 
taking the partial trace with respect to the ancilla variables. 



/ M I -MM* I 

Ut —M* U 2 



Note Physical reasons require that the dimension of J'(M) be (mx k) x (m x k) where 
the ancilla has state space of dimension k. This can always be achieved via direct sum 
with an identity matrix of appropriate size. 

5.3. Mocking Up a Quantum Operation. We now extend our discussion to general 
quantum operations, of which measurement is a special case. We are interested in the 
finite dimensional case. 

Definition 4. A quantum operation is a completely positive map $ : C nxn — > C mxm 
between density matrices that does not increase the trace. 

It is well known that $ must take the form = ^^^EipE* 0. By mocking 
up [9 we mean, similar to the POVM case, coupling the system to an ancilla and find 
an unitary evolution on the combined system such that the reduced state, obtained 
via the partial trace, is Y^i=i EipE*. In other words, given a CP map $, we wish to 
find an unitary dilation. Again Stinespring's theorem ensures the existence of such a 
dilation. But in practice perhaps one would like to obtain a less abstract solution. The 
parametrization of matrix contrations provides one such explicit and easy procedure. 
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Choose the ancilla to be C nm . Any quantum state can be purified and we can 
assume the ancilla is in a pure state, a rank-1 projection, of the special form eoeg, 
where 



1 




e 



The state of the composite system is 



e e ® p 



p 




If U is the proposed unitary evolution, then 



Up U* 



U u U 12 
U21 U22 



p 




TT* TT* 

U ll U 21 

TT* TT* 

u 12 u 22 



U 21 pUt x U 2lP U* 2l 



Tracing out the ancilla, the first system, gives the reduced density matrix Xa™ Ui\pU* v 
Therefore to specify U is to find appropriate operators A, B, and C such that 



U 



is unitary, where T is the contraction 



T A 
B C 



e C 



kxk 



As before, this can be achieved by forming the operator 



J\T) 



I 

Ui 



T D T * 
D T -T* 



I 

U 2 



Notice purfication of mixed states was applied only to the ancilla. There is another 
approach that relies on purification more heavily Namely, one treats the matrix 
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($(Eij))ij, where Eij are standard matrix units, as a state and consider its purifi- 
cation. The differences are, first, in the latter only the range of $ is coupled to an 
ancilla and, second, the lifted map still may not be an unitary evolution. 
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